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^p^l Abstract. Starting from nonlocal symmetries related to Backlund transformation (BT), 

• many interesting results can be obtained. Taking the well known potential KdV (pKdV) 

, equation as an example, a new type of nonlocal symmetry in elegant and compact form which 

comes from BT is presented and used to make researches in the following three subjects: two 
sets of negative pKdV hierarchies and their corresponding bilinear forms are constructed; the 
nonlocal symmetry is localized by introduction of suitable and simple auxiliary dependent 
^ ' variables to generate new solutions from old ones and to consider some novel group invariant 

On , solutions; some other models both in finite dimensions and infinite dimensions are generated 

by comprising the original BT and evolution under new nonlocal symmetry. The finite- 
■ dimensional models are completely integrable in Liouville sense, which are shown equivalent 

to the results given through the nonhnearization method for Lax pair 

o 

(N 



PACS numbers: 02.30.Ik, 11.30.Na, 04.20. Jb 

1. Introduction 

With the development of integrable systems and solion theory, symmetries |[I]-[3l play the 
more and more important role in nonlinear mathematical physics. Thanks to the classical or 
nonclassical Lie group method, Lie point symmetries of a differential system can be obtained, 
from which one can transform given solutions to new ones via finite transformation and 
construct group invariant solutions by similarity reductions. However, little importance is 
attached to the existence and applications of nonlocal symmetries [|2l[3]|. Firstly, seeking for 
nonlocal symmetries in itself is a difficult work to perform. One of our authors (Lou) has 
made some efforts to get infinite many nonlocal symmetries by inverse recursion operators 
ll4l[5]| the conformal invariant form (Schwartz form) (6] and Darboux transformation [lZ1[8]|. 
Moreover, it appears that the nonlocal symmetries are rarely used to construct explicit 
solutions since the finite symmetry transformations and similarity reductions can not be 
directly calculated under the nonlocal symmetries. Naturally, it is necessary to inquire as 
to whether nonlocal symmetries can be transformed to local ones. The introduction of 
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potential [jBJ and pseudopotential type symmetries [|9]-fTTII which possesses close prolongation 
extends the applicability of symmetry methods to obtain solutions of differential equations 
(DEs). In that context, the original given equation(s) can be embedded in some prolonged 
systems. Hence, these nonlocal symmetries with close prolongation are anticipated [il2i - tl4il . 

On the other hand, to find new integrable models is another important application of 
symmetry study. A systematic approach have been developed by Cao [USU TTl to find finite- 
dimensional integrable systems by the nonlinearization of Lax pair under certain constraints 
between potentials and eigenfunctions. Especially in the study of (l+l)-dimensional soliton 
equations, various new kinds of confocal involutive systems are constructed by the approach 
of nonlinearization of eigenvalue problems or constrained flows [fT8l[T9l . It has also been 
pointed that by restricting a symmetry constraint to the Lax pair of soliton equation, one can 
not only obtain the lower dimensional integrable models from higher ones, but also embed 
the lower ones into higher dimensional integrable models [i6l[8l|20l. Here, alternatively, we 
are inspired to act the given nonlocal symmetry on the Backlund transformation (BT) instead 
of Lax pair to generate some other new systems via symmetry constraint method. The related 
work may be adventurous but full of enormous interest. 

In this paper, taking the well known potential KdV equation (pKdV) for a special 
example, we will study the nonlocal symmetry defined by BT. Since the BT reveals a finite 
transformation between two exact solutions of DEs, it must hint some symmetry. For pKdV 
equation, a new class of nonlocal symmetries are derived from its BT, which may give 
more interesting applications than those nonlocal symmetries only including potentials and 
pseudopotentials. The prolongation of the new nonlocal symmetries are found close after 
extending pKdV equation to an auxiliary system with four dependent variables. The finite 
symmetry transformation and similarity reductions are computed to give exact solutions of 
KdV equation. What we want to mention is the process can once lead to two exact solutions 
from one given result due to the Backlund transformation. Moreover, for the pKdV equation, 
some other models both in finite dimensions and infinite dimensions are obtained. The finite- 
dimensional systems obtained here are found equivalent to the results given by Cao [fT6ll . 
which have been verified completely integrable in Liouville sense. This discovery confirms 
that these obtained infinite-dimensional models should have many nice integrable properties, 
which needs our further study. 

The paper is organized as follows. In section II, we present a detailed description 
about the new nonlocal symmetry with BT of pKdV equation. Two kinds of flow equations 
corresponding to the given nonlocal symmetry, i.e. the negative pKdV hierarchies, are 
obtained and their corresponding bilinear forms are also given out. In section III, we extend 
the nonlocal symmetry to be equivalent to a Lie point symmetry of a auxiliary prolonged 
system admitting pKdV equation and its BT. Then the finite symmetry transformation and 
similarity reductions are made to produce exact solutions of pKdV and then KdV equation. 
Section IV is devoted to constructing various integrable systems by means of symmetry 
constraint method. Conclusions and discussions are given in Section V. 
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2.1. BTfor the pKdV equation 

The well-known KdV equation reads 

(Ot + (Oxxx - 6tt)Ct)x = 0, (1) 

where subscripts x and t denote partial differentiation. For convenience to deal analytically 
with a potential function u, introduced by setting £0 = Ux, it follows from equation (1) that u 
would satisfy the equation 

ut + Ujcxx - 3ul = 0, (2) 

which is called potential KdV (pKdV) equation. 

For equation (2), there exists the following BT [[2T| 

Ux + ui,x = -2X+ ^''~p^\ (3) 

Ut + ui^t = 2ux + 2ul^x + 2'^xUi,x-{u-ui){uxx-ui^xx) (4) 

with A being arbitrary parameter. 

Equations (3) and (4) show that if w is a solution of equation (2), so is ui, that is to say, 
they represent a finite symmetry transformation between two exact solutions of equation (2). 

On the other hand, equations (3) and (4) can also be viewed as a nonlinear Lax pair of 
equation (2). For 

mi,a: = -"x-2AH , (5) 

uij = -Ut + 2ul + 2u\j^ + 2uxUi^x-{u-ui){uxx-u\^xx), (6) 

its compatibility condition u\^xt = u\^tx is exactly equation (2). In fact, both equation (3) and 
equation (4) hint that they are all Riccati type equations about uorui, which can be linearized 
by the well known Cole-Hopf transformation 

u = -2—, or Ml = -2^. (7) 
¥ ¥i 

Moreover, by virtue of the dependent variables transformation (|7]), one can convert 

equation (2) into the following bilinear form 

{Dt+DxDr)xi/y/ = 0, (8) 
meanwhile it leads equations (3) and (4) to 

(d2-A)v/-V^i=0, (9) 
{D,+Dl + 3?iDx)xir-x{ri=0, (10) 

where the Hirota's bilinear operator D"^Df is defined by 



51-57) [jrir) ''(^•'>*(^-'' 



x'=x,t'=t 
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2.2. The nonlocal symmetry from Backlund transformation 

For equation Q with its BT ([3]) and dH), considering the invariant property under 

we may find substantial possible nonlocal symmetries and a special case is presented and 
studied as follows. 

Proposition 1. The pKdV equation (2) has a new type of nonlocal symmetry given by 

a = exp{ J u — uidx) , (11) 

where u and ui satisfy BT (3) and (4). That means a given by (fTTI) satisfies the following 
symmetry equation 

Ot + Oxxx - 6uxOx = 0. (12) 

Proof: By direct calculation. 

On the other hand, we let the bilinear pKdV equation ([8]) be invariant under the 
transformation i// — )• i// + eay/, which produces the corresponding symmetry equation 

{Dt+DxDt)o^-xir = 0. (13) 

The Cole-Hopf transformation u = —2^ between equation (2) and its bilinear equation ([8]) 
determines a symmetry transformation for a and 0^, saying 

^ ^ 2v/,gy _ 2gy,^ 

Taking equations (fTTI) and ([7]) into equation (fT4l) . we obtain a class of nonlocal symmetry for 
equation ([8]) 



Correspondingly, it gives the following proposition for equation ([8]). 

Proposition 2. The bilinear pKdV equation ([8]) has the nonlocal symmetry expressed by 
(fTSl ), where i// and i/Zi satisfy bilinear BT © and (fTOl) . 

Proof: One can directly check that a^, given by (fTSi) satisfies symmetry equation (fT3l) under 
the consideration ([91) and (fTOl) . 

2.i. Two sets of negative pKdV hierarchies 

The existence of infinitely many symmetries leads to the the existence of integrable hierarchies 
and with the help of infinitely many nonlocal symmetries, one can extend the original system 
to its negative hierarchies [|22ll23il . Here, starting from the nonlocal symmetry ([TT]) related to 
BT of equation (2), we would like to present two sets of negative pKdV hierarchies and their 
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corresponding bilinear forms are also constructed only by the transformation (|7]). 
Case 1. The first kind of negative pKdV hierarchy can be obtained, reading 

N r. 

"f-w = -I^exp( / a-a,djc), (16) 
(=1 •' 

u, + u,,, = -2h+^^^-^, /=1,2,...,A^, (17) 

where A, is arbitrary constant. 

In particular, when = 1 , one has the first equation of negative pKdV hierarchy, namely 

luxxtUt-'^UxU^ -u^^f-AXiu^ = Q. (18) 

Here we have instead t-\ with t for simplicity. It is well known that the first negative flow in 
the KdV hierarchy is linked to the Camassa-Holm equation via a hodograph transformation 
ll24l or can be reduced to the sinh-Gordon/sine-Gordon/Liouville equations [|25l . Here we 
transform equation (fTSi) into sine-Gordon and Liouville equations. 

In fact, by setting = fi{x,t) = —Uf, we can rewrite equation (fTSl) in the form 

^=(-k + Mj)^ (19) 



2/3 4/3V; 
which can be integrated once with respect to x to give 

/3(ln/3),, + i82 = /3o(0, (20) 
where /3o(?) is an arbitrary function of t. 



As it is reported in Ref. [|24l . for non-zero Po{t), one can rescale j8 to ^ /3o(?)/3, redefine 
tast/ a//3o(0 and set /3 = exp(/T]) to give the sine-Gordon equation 

rixt = smri, (21) 

while for j8o(0 = 0, by setting /3 = — exp(77), equation (|20|) becomes the Liouville equation 

T?xr = er (22) 

Remark 1. The quantity — ^ + fp: = ^ in the right hand side of equation (fT9l) can be given 
in terms of a Miura transformation 

A = -e,-G\ e = ^. (23) 

Furthermore, by virtue of the dependent variable transformation 

u = -2—, Ui = -2^, {i=l,2,...,N) (24) 
the negative pKdV hierarchy (fT6l)- (fT7l) is directly transformed into its bilinear form 

N 

D,Dr^yfxif=Y,¥f. (25) 

(=1 

{Dl-?ii)xif-xifi = 0, /=1,2,...,M (26) 
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Case 2. For the nonlocal symmetry (fTTl) being dependent with parameter A , we may derive 
the second kind of negative pKdV hierarchy by expanding the dependent variable in power 
series of A . In this case, we have 

1 f exp{J u — uidx) 



Ut_ 



Nl 



[U — U[ 



A=0 



(27) 



(28) 



Under the transformation u = — ^ and ui = — the negative pKdV hierarchy (|271) -(|28]) 



becomes 



1=0 



(29) 
(30) 

(31) 



1 fd^\iff\ 

(d2-A)v/-V^i=0. 
Let v/i = V^i (A) have a formal series form 

oo 

;=0 

where y/, is A independent. Then (|29l)- (|30l) can be rewritten as 

A' 

D,Dt^ i/A ■ i/A = £ fk¥N~k, (32) 
DlY-Vk = VVk-i {k = 0,l,...,N) (33) 

with =0. 

The negative pKdV hierarchy in bilinear form (I32l)-(l33l) is just the special situation of 
the bilinear negative KP hierarchy for dy — O 'm Ref. L22J. From this observation, we have the 
following remark: 

Remark 2. The second negative pKdV hierarchy shown by (|27l) -(l28l) is a potential form of a 
known negative KdV hierarchy given by other methods, say, the inverse recursion operator dU, 
Lax operator [|23l . and the Guthrie's approach |[26ll . 



3. Localization of the nonlocal symmetries 

We know that the Lie point symmetries |[2l[3l| can be applied to construct finite symmetry 
transformation and group invariant solutions for DEs, whereas the calculations are invalid for 
the nonlocal symmetries. So it is anticipant to turn the nonlocal symmetries into local ones, 
especially into Lie point symmetries. In order to make the nonlocal symmetry localized, one 
may extend the original system to a closed prolonged system by introducing some additional 
dependable variables [[T2] - [T4| to eliminate integration and differentiation. 

Fortunately, starting from the nonlocal symmetry (fTTl) . the prolongation is found to be 
closed when another two dependent variables v = v(x, t) and g = g{x, t) are introduced by 



Vx = u-ui, 
gx = e", 



Vt = 2{u-ui){ux-2X) -2mxx, 
gt = -e^'[2u^ + SX - (m-mi)^]. 



(34) 
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Now the prolonged equations (2), (3), (4) and (|34l) contain four dependable variables u. 
Ml, V and g, whose corresponding symmetries are 

1 



Ou = e", a„, = 0, Oy = g, Og = -g . (35) 



gt = {g;x}g, + 6Xg,., {g;x}^^--^, (36) 



.2 

Remark 3. What is more interesting here is that the symmetry shown in (|35l) implies the 
auxiliary dependent variable g satisfies 

which is just the Schwartz form of the KdV (SKdV) equation (1). This may provide us with a 
new way to seek for the Schwartz forms of DEs, especially for the discrete integrable models, 
without using Painleve analysis. 

Due to (l35l) . the symmetry vector of the prolonged system has the form 

y = eV+0^+^^ + -g2 (37) 
ou oui ov 2 og 

Then, by solving the following initial value problem 

du ,7 dU\ dv _ dg 1 _2 

de~^ ' 1e~ ' d£~^' de ~ 2^ ' (38) 

u\e^O = U, Mi|e=o = "l, v|e=o = V, ^|e=0=,?, 
the finite transformation can be written out as follows 

2e V 2 2 

u = u + - e, ui=uu v = v + 21n- , g = - g. (39) 

2-£g 2 -eg 2 -eg 

Remark 4. The original BT (3) and (4) in itself suggests a finite transformation from 

one solution u to another one u\ and then the new BT (|39l ) obtained via (fTTI) will arrive 

at a third solution U. Actually, the finite transformation (l39l) is just the so-called Levi 

transformation If27l . The result of this paper shows the fact that two kinds of BT possess 

the same infinitesimal form (fTTI) . 

Now by force of the finite transformation (|39l ). one can get new solution from any initial 

solution. For example, it is easy to solve an initial solution of prolonged equation system (2), 

(3), (4) and dM]), namely 

M = c, Ml = c-l-2 V^tanh^, v = — ln(tanh^ ^ — 1), 

sinh(2C) X . nr, . . (^^^ 

4 VA 2 

Where A, c and cq are three arbitrary constants. 

Starting from this original solution (|40l) . a new solution of equation (2) can be presented 
immediately from (l39l) : 

8 VXe cosh^ C, ^ ^ . ^ 

u = c = = , (41) 

8 yi- £[sinh(2C) - 2 - \2Xt - 2co)] 

which then gives the corresponding solution of KdV equation 

_ ... [cosh(2g) + 1]£ + v/Isinh(2C) \A^e{x-\2Xt- 2co)] 
[8 VI-e(sinh(2C)-2 12A?-2co))]2 
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with C = y/X{~x + 4Xt). 

Besides obtaining new solutions from old ones, symmetries can be applied to get special 
solutions that are invariant under the symmetry transformations by reducing dimensions of 
a partial differential equation. To find more similarity reductions of equation (2), we will 
study Lie point symmetries of the whole prolonged equation system instead of the single 
equation (2). Suppose equations (2), (3), (4) and (l34l) be invariant under the infinitesimal 
transformations 

M— T-M + ea, Ml — )► Ml + ecTi, v^v + £02, g^g + eo^, 

with 

=X{x,t,u,ui,v,g)ux + T{x,t,u,ui,v,g)ut-U{x,t,u,ui,v,g), 

01 =X{x,t,u,uuv,g)ui^^ + T{x,t,u,ui,v,g)ui^t-Ui{x,t,u,uuv,g), ^^^^ 

02=X{x,t,u,ui,v,g)vx + T{x,t,u,uuv,g)vt-V{x,t,u,ui,v,g), 

03=X{x,t,u,uuv,g)gx + T{x,t,u,u\,v,g)g,-G{x,t,u,ui,v,g). 

Then substituting the expressions (|43l) into the symmetry equations of equations (2), (3), (4) 
anddMl) 

Ot + cTxxx - 6m^o:^ = 0, 

0\^x + Ox - [o - 0\){u- U\) =0, 

0\t - Ojcxx + 2{u -U\)Oxx + 2(C7 - Ox)Uxx - [4A + (m - Uif - 2Ux]Ox 

+2{o-Oi){u-ui){2X-Ux) =0, 

0-2,x - CT + (7i = 0, 

02t + 20xx + 2{ui-u)Ox + 2{oi -o){ux-2X) = 0, 
cy3,x-e''a2 =0, 

1 , 

03t + 2e''[Ox + {Ui - u){0 - Oi) - -{u - U1Y02 + i^X + Ux)02\ =0, 

and collecting together the coefficients of partial derivatives of dependent variables, it yields a 
system of overdetermined linear equations for the infinitesimals X, T, U ,U\,V and G, which 
can be solved by virtue of Maple to give 

X{x,t ,u,u\,v,g) = c\{x+ \2Xt) + C5, 
T{x,t,u,ui,v,g) = 3cit + C2, 

U{x,t,u,ui,v,g) = —ci{2Xx + u) +2c4e^ + C3, ^^^^ 

Ul{x,t,U,Ul,V,g) = -Cl{2Xx + Ul) +C3, 

V{xj,u,ui,v,g) = -ci +2c4g + ce, 

G{x,t,U,UuV,g) =C4^2 + C6^ + C7, 

where c,(z = 1...7) are seven arbitrary constants. When ci = C2 = C3 = C5 = = C7 = 0, the 
reduced symmetry is just (l35l) . 

To give the group invariant solutions, we would like to solve symmetry constraint 
conditions a = and a, = 0(z = 1,2,3) defined by (l43l) with (l45l) . which is equivalent to 



(44) 
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solve the following characteristic equation 

dx dt du 



ci{x+12Xt) +C5 3cit + C2 —ci{2Xx + u) + 2c4e^' + C3 
dui dv dg 

-Ci(2AjC + Mi) +C3 -Cl+2c4g + C(, C4g^ + C(,g + C']' 

Two nontrivial similar reductions under consideration C4 are presented and substantial 
group invariant solutions are found in the follows. 
Case 1: ci and — 4c4C7 0. 

Without loss of generality, we let ci = 1 . For simplicity, we introduce arbitrary constants 
a4 and a-j to replace C4 and c-j by al = c-^ — 4C4C7 and a-j = — a|/(16c4), then after solving 
equation (|46|) . we have 

U = -Xx + 3X^t + C3+C5?i-3c2X^ + i3t + C2y^U{^) 

'''■exp(y(^)-G(^))tanh5], 



4a-j 



ui = —Xx + 3X ? + C3 +C5A — 3C2A + 



{3t + C2)V (47) 



v = -iln(3r + C2)-G(^)+y(^)-21ncosh5, 

g = [tanh5H J 

a4 a4 

with B = a4i3G{^)+ ln(3ci? + C2) )/6 and ^ = (jc - + C5 - 6C2A) /(3? + C2) ^ . 

Here, U{^), t/i(<^), V{^), G{^) and represent five group invariants and substituting 
(1471 ) into the prolonged equations system gives the following reduced equations 

1 a^ 

H,, = ^+4ajH^-tH (48) 

2H ^ 32a^H' 



f/(0 = f/i(^)-^, y(^) = G(^)-ln(//), G^(^; ^ 



(49) 



4^7^ 

with H = H{^). One can see that whence H is solved from equation (l48l) . two new group 
invariant solutions w and ui of equation (2) would be immediately obtained through equations 
dH]) and (Ell). 

Moreover, by making a further transformation [f28l 

^(^) = ^(^§+/'' + |), ^ = (50) 

equation (l48l) can be converted into the second Painleve equation Pn, reading 

P^^=2P^ + ^P + a, (51) 

with a = — (a4 + l)/2. Now, every known solution of Pn (ISTI) will generate two new group 
invariant solutions of equation (2), and then two new solutions of KdV equation (1) denoted 
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as (Oi and (O2 can be given directly after one derivative with respect to x for ui and u 

coi = ^ ^{P^ A, (52) 

1 , a| 9 ,2aAP a\ ^ laAP 7, 

£02 = ^[--^scch^Ri + (-^ - t4) tanhi?i + - p2 + A, (53) 

(3? + C2)^ F F ^ 

where 

F^F(^)=2/'^+2/'2 + ^, 7?i = ia4[ln(3? + C2) + 3G(^)], G^(^) = ^^^i^^^, 

and P satisfies Pn dlB with a = -(a4 + l)/2. 

It is known that the generic solutions of Pn are meromorphic functions and more 
information about Pn is provided in Ref. [29|, saying: (1) For every a = N E Z, there exists 
a unique rational solution of Pn; (2) For every a = N + j, with N E Z, there exists a unique 
one-parameter family of classical solutions which are expressible in terms of Airy functions; 
(3) For all other values of a, the solution of Pn is transcendental. 

For example, when a = 1 (a4 = —3), Pn (|5TI) possesses a simple rational solution 
P{£,) = -l/^, which leads the solutions ^ and ^ to 

2 

[x — oAt + C5 — 0C2Ay 

and 

m2 = -[x^- 36tx^ + (540^2 - 6)/ - (4320? ^ - 1 68? - 2)^3 + 36? (540?^ - 48? - 1 )x^ 
-(46656?^ - 7776?^ - 216?^ - 144? - 12)jc + 46656?^ - 12960?^ - 432?^ 
-720?2 - 48? + l]/[x^ - ISx^t + 108jc?2 - (6? + 1) (36?^ - 6? - 1)]^ (55) 

In the formulation (|55l) . we have made C2 = 0, C5 = and A = 1 because the original 
expression is much too complicated. The simple rational solutions of PII will yield abundant 

rational solutions of KdV equation. 

1 
2 



When a = k {a4 = —2), Pn (l5TI) has a solution expressed by Airy function 



_.3Ai(l,-24^)-v/3Bi(l,-24^)_ (56) 
3Ai(-2-3^)- V3Bi(-2-3^) 
For simplicity, we convert equation (l56l) into the equivalent form 

where J(n,(^) is the first kind of Bessel function. Substituting (l57l) into (l52l) and (|53l) with 
C2 = C5 = and A = 1 (or else the formulae are too long to written down here), two exact 
solutions of KdV equation are obtained as follows: 

, ;c-12? Jc-6?!, 

»I =— if- <^*' 

4 = -f/a (59) 



Nonlocal symmetries related to Bdcklund transformation and their applications 



11 



with 



^ = 32 Vt{x - 6t)\x{x -6tf- \2t)'S\+x{x - 6tfil] + 128 \/6t{x -6t)^ ■ {j\h 
+J^Ji ) + 96 Vt{x - 6t)\x{x -6tf- 4?) Ji J2 + - etft2^\] -11-1^ ^t 
■{x - ^tf ■ {{x{x - 6tf - At) ■ fi+x{x - 60^J2] + 256 V6t{x - 6?)^ Ji J2 
~\92V2>-2h{x- 6t)i ■ {j]h + J^Ji) +54-2^ -x V~t{x - 60^(J? + J2) 
-144 • 2h Vt{x - + 144 V3 ■ 2s^(x - 6f)3Ji J2 - 21x ^t, 

and 

Q. = 6f5 [23 {x - 6tf{i\ + J^) - 3] [8 ■ 23 (x - 6tf{)\ - J^)^ - 12 ■ 2? (jc - 60^(J? + J2) + 9] , 

where we denote Ji = J(l, ^^^), J2 = J(-f , ^^^). 

Then continue to do the same, sequences of rational solutions and Bessel (Airy) function 
solutions for KdV equation will be easily constructed. Furthermore, by selecting suitable 
parameters in this kind of similar reduction, we may discover more unknown exact solutions 
among interaction solitons and Painleve waves of KdV equation. 
Case 2: ci = and C2 0. 

Firstly, it is convenient to replace C4 and C5 with a^ and k hy a\ = c^ — Ac^c-] and 
k = C5/C2, and it follows the results from equation (l46l) . saying 

u = % + U{z) + ^G(z) - Mzil eV(z) tanh[^l(i±^] , 

C2 C2 a4C7 2C2 

Ml = —t + UAz), 

2C7 , ^M^{t + G{z)),, ^^^^ 

g = [C6 + a4 tanh( ^^)] , 

a^ ZC2 

v = y(z)-21ncosh[^^l^^±^] 

2C2 

with z = X — kt. Substituting (|60l) into equations (2), (3), (4) and (|34l) and redefining the 
parameters for the sake of simplicity, we notice that the new group invariants U{z), Ui{z), 
V{z) and G{z) are subject to 

W^-alW^^-a^W^+asW^-ajW = 0, (61) 



2^^ 4^ 



y(z) =ln(W^), 

G,{z) = - 
a-j 

With W = Wiz) and a2 = — — as = ^ — f as = 2k- 12A, aj = 

\ ' ^ 04^7 046-7 
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Remark 5. The case ci = here is interesting. From equation (|6T| ). we know that W can 
be expressed as an elliptic integration and can be expressed by means of Jaccobi elliptic 
functions. Whence W is fixed from (|6TI) . all the other quantities are given simply given by 
differentiation or integration. The first equation of (l60l) implies the important byproduct, the 
explicit exact interaction between cnoidal periodic wave and kink soliton. 

A simple example of this case can be obtained by using the simplest Jacobi Elliptic 
function expansion method which leads to 

W{z) = ^[sn{^,n)-l] (63) 
4aj Aain 

with the constraint conditions a^ = "^^5,^2^" ^^'^ = '^tirfij'^ equation (|6TI) . where n is 
the modulus of the Jacobian elliptic function sn. 

After solving equation (|62l) with the given solution (|63l) and taking the results into (|60l) . 
two exact solutions of the KdV equation are obtained 

al(l-n^) ai(5n2-l) ^ 
(03 = — -r-r: — ^+ ^ . ' -A, (64) 
16a^n2y + l 64aln^ 



(04 = ^^ ^tanh2i?2 + ^-^^ tanhi?2 + - 



32a2 16na2 a^ 

al{n^+l] 



n^a2 



-64A, (65) 



7 2 ' 



where we have 



a\{j?-\){t^a(,) I f In 

K2 = ^ o ■ r T 



2 

64a3n2 ' A"\2n^Y^ + 2nVY^^^fi^Y^^-n^ 
1 1 /a3(192AaV-5aV + a2)/' agx 



2 Jo v^r^Vl-n2f2 V 128n3a3 4a2n 

and (32, as, and A are four arbitrary constants. 



4. Integrable models from nonlocal symmetry with Backlund transfer- mation 

To find new integrable models is another important application of the symmetry study. 
Symmetry constraint method is one of the most powerful tools to give out new integrable 
models from known ones. Especially, casting symmetry constraint condition to Lax pair of 
soliton equations, one can obtain many other integrable models. In this section, we would like 
to combine the nonlocal symmetry with BT of pKdV equation to give some integrable models 
both in lower and higher dimensions. 

Let every pair {u,Ui) (i = l,2,...,N) satisfy the following BT 

Ux + M/,x = -2 A, H , (66) 

ut + Ui^t = 2ul + 2u\^ + 2mxM/,x - (m - Ui) [ujcx - Ui^xx) , (67) 
and the corresponding nonlocal symmetry of u reads o' = exp(/ u — Ujdx) for / = 1,2, A^. 
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4.1. Finite-dimensional integrable systems 

In general, every one symmetry of a higher dimensional model can lead the original one to its 
lower form. Now, considering 

A' 

(68) 



It n 

= 52 a,exp( u — Uidx) 
i=i 



as a generalized symmetry constraint condition and acting it on the x-pait of the BT (3), we 
firstly give the finite dimensional (A^+ 1) -component integro-differential system 

N . 

Ux = ^ ai exp( I u — Uidix) , 

'=1 (69) 

+ = -2A,:H , 2 = 1,2,...,A^, 

where every a, and A,- are arbitrary constants. For further simplification, making m,- = 
u — {\nwix)x, then the constraint condition (|68l) becomes 

N 

M = 52 ^mWm, (70) 

which transforms (l69l) into the N-component differential system 

N 

^WijcxxWix - 4( 52 amWmx)wjx " " "^^i^fx = 0, « = 1, 2, A^. (71) 

Taking Si = Wix, we rewrite equation (TtTI) as 

2*w;c*i -4(5^ ^ilm^m)*? " sjx - ^^isj = 0, / = 1 , 2, . . . , A^. (72) 

m=l 

,2 



Making Si = bigf, equation 0721) is equivalent to the downward integrable system 

N 

qixx - ( 5^ Cmqlr)qi - kqi = 0, / = 1 , 2, . . . , (73) 

with c, = a/Z?, being arbitrary constant. 

On the other hand, by the same symmetry constraint (168]) and the ?-part of the BT (4), we 
can construct another set of integrable system 

A' 



(74) 

Mf + M,- , = lul + 2m^,^ + 2UxUi^x -{U- Ui) [Uxx - Uf^xx) , z = 1 , 2, . . . , A^. 

Considering the similarity transformations of dependent variables done in x-part, after a series 
of tedious substitutions, equation (1741) becomes 

N N N 

qiqixt - quqix + (2 52 ^^qli - ^^dqlx - ^qtqix 52 '^mqmqmx -'^qKY^ Cmqlrf' 

m=l m=l m=l 

N N 

+2Xiqj 5^ c,nql, + 4X1 qj + 2'?^ 51 c„ (^^_, + qmqmxx) =0, / = 1 , 2, . . . , A^. (75) 

ni=l m=l 
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Taking equation (1731 ) into account, equation (1751) can be integrated once about x to give A'^- 
component integrable system 

N 

qu = -2 ^ Cmqmqmxqi + 2 Cmqliqix-'^^iqix, i = 1,2, ...,N. (76) 

m=l m=l 

In fact, equations (1731) and (1761) are essentially the canonical equation (Fq) and (Fi) 
respectively [16], saying 

A' 

(Fo) : qix = Pi, Pix = ( (^inqm)qi + hqi- (77) 

m=l 

qu = -2{Zm=l CmPmqm)qi + 2{Xm=\ Cmqm)Pi " 4A/P,-, 
iFl):{ Pit =2{Y!i=lCmPmqm)Pi-2{Y.'l=\^mPm)qi-^^fqi (78) 
-2Hl!^n=\ Cmqi)qi - 2{Lnr=l ^mCmq^qi- 

It should be stressed here that the finite integrable systems (1731) and (1761) reobtained via 
this way are just the remarkable results given by Cao in Ref. [ [T6l through the nonlinearization 
method, both of which have been proved completely integrable in Liouville sense. Thanks 
to these finite integrable systems ([771) and (1781) . the original high dimensional KdV equation 
would be solved. 



4.2. Infinite-dimensional integrable systems 

For getting some higher dimensional integrable models, one may introduce some internal 
parameters |l6l|8l|20l|. Here, we would like to use the internal parameter dependent symmetry 
constraints on BT to construct two sets of infinite-dimensional integrable systems. 

It is obvious that equation (2) is invariant under the internal parameter translation, say y 
translation, so we can view 

N r. 

Uy = '^ ai exp( u — Uidx) (79) 
i=i 

as a new symmetry constraint condition. 

Firstly, imposing (1791) on the x-part of the BT (3) yields a (l+l)-dimensional (A'^+ 1)- 
component integro-differential system 



N . 

My = ^ ai exp( u — Uidx) , 
i=i 



iu — Ui)^ 

Ux + Ui^x = -lAi^ , i^l,2,...,N, 



(80) 



where A,- and ai {i = l,2,...,N) are constants. 

By the transformation Ui = u — {\n(^iy)x, equation (1791) becomes 

N 

u = Y^ai^i, (81) 



i=l 
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which then converts (|80l) into the following (l+l)-dimensional N-component differential 
system 

2^ixxy4iy - 4 ( £ a,n(^,n, ) 0j - 0,2^,, - 4A,(^)2 = 0, / = 1 , 2, . . . , A^. (82) 

Alternatively, combining the constraint condition (1791) with the ?-part of the BT (4) will 
produce a (l+2)-dimensional system about x, y and t, reading 

A' 



f 

Uy = £ ai exp( / M — M,djc) , 



(83) 

u, + J = 2m^ + 2u\^ + 2m;,m,,x - (m - Ui) (mxx - M/,xx) , «■ = 1 , 2, . . . , A^. 
Using the same transformation and equation (|82l) . equation ([83] ) is transformed into the 
following A^-component system 

A' A' 



^ixyt^iy - ^ixy^iyt " 2( a,n^nixx)^iy^ixy + ( Cl,n^mx " 2A,-)^,^ 



m=l m=l 
N N N 

+2[( a,n^rnxf + 2A/ - + 4Af ] = 0. (84) 

m=l m=l m=\ 

It should be noted that the integrability of the infinite-dimensional systems (f82l) and 
(f84l) obtained in this way is not quite clear. The finite-dimensional models obtained here are 
completely integrable, that strongly suggests these infinite-dimensional models should have 
many nice integrable properties. It will be of much interest to investigate the integrability of 
these models in the further work. 



5. Conclusion and discussions 



In this paper, we have shown that combining nonlocal symmetries with BTs can result in 
many diverse applications. The main new progresses made in this paper in the general aspect 
of integrable systems are: 

(i) . The BTs are used to find nonlocal symmetries; 

(ii) . Different types of BTs may possess same infinitesimal forms and then new types of BTs 
may be obtained from old ones; 

(iii) . New integrable (negative) hierarchies can be obtained from nonlocal symmetries related 
to BTs; 

(iv) . New finite dimensional integrable systems can be obtained from BTs and related 
symmetry constraints and reductions. And then the original high dimensional model can 
be solved from lower dimensional ones because of the existence of nonlocal symmetries 
depending on BTs ; 

(v) . The exact interaction solutions among soli tons and other complicated waves including 
periodic cnoidal waves and Painleve waves are revealed which have not yet found for any 
integrable models because it is difficult to solve the original BT (or Darboux transformation) 
problem if the original seed solutions are taken as the cnoidal or Painleve waves; 

(vi) . The localization procedure results in a new way to find Schwartz form of the original 
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model which is obtained usually via Painleve analysis for the continuous integrable systems. 
The method may provide a potential method to transform a discrete integrable systems to 
Schwartz forms because usually the BTs of discrete integrable models are known. 

The above progresses are realized especially for potential KdV (pKdV) equation. For 
pKdV equation, it possesses a new class of nonlocal symmetry resulting from its BT. Since 
this BT is of Riccati type, more information about its bilinear forms is learned via the Cole- 
Hopf transformation. Based on the new nonlocal symmetry with internal parameters, we 
construct two sets of negative pKdV hierarchies and fulfill their corresponding bilinear forms. 

In order to extend applicability of nonlocal symmetry to obtain explicit solutions of 
KdV equation, we introduce another two auxiliary variables v and g to form a prolonged 
system with u and wi, so that the original nonlocal symmetry can be transformed to a Lie 
point symmetry of the new equations system. Then what follows naturally are Lie-Backlund 
transformation and two kinds of novel similarity reductions. By virtue of two kinds of BTs, 
the solitary wave solutions of KdV equation are obtained through the transformations of the 
trivial solutions. Concerning the complete Lie point symmetries of the prolonged system, we 
achieve rich group invariant solutions including rational solution hierarchy, Bessel function 
solution hierarchy and periodic function solutions. 

The nonlocal symmetry has also been devoted to construct various new integrable 
systems by symmetry constraint method. Applying nonlocal symmetry on the BT of pKdV 
equation, finite-dimensional integrable systems are given, which are found equivalent to the 
excellent work done by Cao [16]. Moreover, the introduction of an internal parameter as new 
independent variable helps us to build two sets of infinite-dimensional models. 

We believe that both the negative pKdV hierarchies and two sets of infinite-dimensional 
models obtained in the paper should have many nice integrable properties. For the completely 
integrable finite-dimensional models, one may consider their algebraic geometry solutions to 
achieve related solutions of KdV equation. For the localization of nonlocal symmetries, it still 
remains unclear what kind of nonlocal symmetries must have close prolongations and can be 
applied to construct exact solutions. It is quite reasonable that these matters merit our further 
study. 
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